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The working principle of magnetic force microscopy and scanning SQUID microscopy is introducing a mag- 
netic source near a superconductor and measuring the magnetic field distribution near the superconductor, from 
which one can obtain the penetration depth. We investigate the magnetic field distribution near the surface of a 
magnetic superconductor when a magnetic source is placed close to the superconductor, which can be used to 
extract both the penetration depth A L and magnetic susceptibility x by magnetic force microscopy or scanning 
SQUID microscopy. When the magnetic moments are parallel to the surface, one extracts A L j -\/l - Any- When 
the moments are perpendicular to the surface, one obtains A L . By changing the orientation of the crystal, one 
thus is able to extract both^ and A L . 
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I. INTRODUCTION 

Superconductivity is well characterized by two length 
scales. The coherence length £ describes the rigidity of the 
phase coherence and the penetration depth Al characterizes 
the response to electromagnetic fields. The penetration depth 
is directly connected to the superfluid density and the pairing 
symmetry, thus its measurement is crucial for the understand- 
ing of new discovered superconductors. The coherence length 
can be measured from the upper critical field H& for type II 
superconductors. There are many well-developed method to 
measure the penetration depth 1 , such as the magnetic force 
microscopy (MFM) and scanning SQUID microscopy (SSM), 
which are the main focus of the present study. 

In MFM and SSM, a magnetic source is placed near a 
superconductoi'^H. \ n MFM, the source magnetic field is gen- 
erated by a small magnetic tip, which can be modeled as a 
point dipole. In SSM, the source field is generated by a cur- 
rent loop. Due to the exclusion of magnetic field by the su- 
perconductor, the magnetic field outside the superconductor 
is modified compared to that without the superconductor. The 
exclusion thus causes repulsion between the magnetic source 
and superconductor. In MFM, the resulting magnetic field dis- 
tribution is measured by the force between the MFM tip and 
superconductor. In SSM, the magnetic field is measured by 
a SQUID. From the measured magnetic field, one can extract 
the penetration depth by fitting to theoretical expressions. For 
non-magnetic superconductors, the magnetic field distribution 
was calculated in Refs.|5K9]for isotropic superconductors and 
in Ref. 10 for anistropic superconductors. 

Recently there is growing interest to apply both MFM and 
SSM to magnetic superconductors, where magnetic ordering 
coexists with superconductivity^E] When a magnetic field 
induced by a source is applied to the magnetic superconduc- 
tors, it polarizes the magnetic moments near the surface of the 
superconductors, which gives additional contribution to the 
magnetic field outside the superconductors. The polarization 
is characterized by the magnetic susceptibility x m the linear 
response approximation. For instance, in MFM, the polariza- 
tion lowers the energy of the whole system, thus gives attrac- 



tion contribution between the MFM tip and superconductor in 
additional to repulsion due to the screening of magnetic field 
by superconductors. The magnetic field distribution outside 
the magnetic superconductor thus depends on x an d ^l- It is 
still an open question what information can be extracted by 
MFM and SSM in the case of magnetic superconductors. Re- 
cently Kirtley et. al. studied the SSM response in isotropic 
paramagnetic superconductors. 14 The effects of the isotropic 
paramagnet in this case are two folds. First it reduces the pen- 
etration depth according to Al ^1 - 4nx- Second, it changes 
the boundary condition. In isotropic paramagnetic supercon- 
ductors, the magnetic field outside depends on Ail -y/1 - 47ry, 
and one cannot extract both Ai and x from SSM measure- 
ments. 

The magnetic superconductors usually have anisotropy in 
magnetic structure. The polarization depends on the orien- 
tation of the magnetic source with respect to the anisotropy 
of the magnetic structure. By changing the orientation of the 
crystal, it is possible to obtain both Ai and x- 

In this work, we investigate the magnetic response in mag- 
netic superconductors based on the London approach both in 
the Meissner state and mixed state. For the magnetic moments 
parallel to the surface of superconductor, the magnetic field 
outside the superconductor depends on A L / ^ 1 - 4-nx when 
the separation between the magnetic source and superconduc- 
tor is much larger than Al. For the moments perpendicular to 
the surface, it depends only on Al. By changing the orienta- 
tion of the crystal, one thus can obtain both the bare penetra- 
tion depth and the magnetic susceptibility. 



II. MODEL 

In this section, we derive the magnetic field distribution in- 
side and outside the magnetic superconductor. A schematic 
view of the setup is shown in Fig. [T] To be specific, we 
consider the case with an easy-axis anisotropy in magnetic 
structure, which is most commonly encountered in magnetic 
superconductors 11 13 . The penetration depth is assumed to be 
isotropic. We consider two cases with magnetic moment par- 



2 



(a) 


z ' 




Vacuum 


\ 


source 

r 


, _ n 


a 




Magnetic 


moments 



(b) 



Vacuum 



0, 



77777777. 

t t t t 

Magnetic 



777777777* 

t t t t 

moments 



Superconductor 



Superconductor 



FIG. 1. (color online) Schematic view of a magnetic source placed 
on top of a magnetic superconductor. The magnetic subsystem has 
easy-axis anisotropy. The easy axis is parallel to the surface in (a) or 
perpendicular to the surface in (b). 



allel to the surface Fig. [TJa) and perpendicular to the surface 
Fig. [3b). 

The magnetic field outside the superconductor is given by 



V x (B - 4ttM s ) = 0, 

where M s is the magnetization in the source. Using V • B 
we can rewrite Eq. ([TJ as 



V 2 B 



-47rVx VxM, 



(1) 
= 

(2) 



Inside the superconductor, we use the London approxima- 
tion which is valid when the penetration length is much 
larger than the coher ence length as realized in most magnetic 
superconductorjnSEl 



A]y x V x (B - 4ttM) + B = <5>a5(x)5(y)i, 



(3) 



where z is the unit vector along the z axis and (Do = hc/(2e) is 
the flux quantum. We assume the vortex density is small when 
the applied magnetic field is much smaller than H c2 , and we 
only consider a single vortex at (x,y) = (0, 0). 

The magnetic field outside and inside the superconductor is 
connected through the boundary conditions at the interface z = 
0. The normal component of B is continuous at the interface 



B z (z = 0") = B z (z = + ). 



(4) 



As the induced surface supercurrent is finite at the interface, 
from the Maxwell equations we have the boundary condition 
for the tangential component B xy 



BAz = 0-) 
B y (z = 0-) 



•4tt 



MAz = 0-) 
M y (z = Or) 



By(Z = + ) 1 • 



The general solution to Eq. Q can be written as B = B t + 
B2, where Bi and B2 are solutions of Eq. Q with M, and 
without M s respectively. 82(1") can be written as 



B 2 (r) 



(2tt) 3 / 2 



dk x dk y B 2 (k) exp [ik^ ■ r 2d - k 2d z] , (6) 



where k^ = (k x ,k y ) and = (x,y). To solve the London 
equation Eq. ([3]), we need to know the magnetic structure of 
the magnetic superconductor. In the following subsections, 
we treat the cases with magnetic moments perpendicular and 
parallel to the interface separately using the linear response 
approximation for the magnetization. 



A. Magnetic moments parallel to the interface 

We assume the easy axis is along the x direction. The mag- 
netization can be written as M x = xB x , which is valid when 
M x <k Mo with Mo being the saturation magnetization. The 
solution to Eq. ([3]) can be written as B = B3 + B4 with B3 ac- 
counting for the magnetic fields induced by the vortex and B4 
being the solution to Eq. ([3]) in the absence of vortices. Since 
the vortex is along the z axis, we have B^ x = #3,y = and 



On 



B 3 Jk x ,k y ) = ■ 

2^4 + 1) 

B4 in the Fourier space is given by 

(k 2 ^ + 1) Bi>x - AnX z L [k 2 y + k 2 z ) x B 4 , x = 0, 



(k 2 4 + l)( ) + 4nAlk x [%\ A B 4 >=(). 



From Eq. ([8]), we obtain 



Anxk 1 - A L 2 - k 2 



2:1 



l-4nx 



From Eq. |9]), we obtain 

k x ky 



B4,y 



kl+k 1 



B\ x , and B\ z 



kl + k} 



B4,; 



(7) 



(8) 



(9) 



(10) 



(11) 



Using the boundary condition Eqs. Q and we have for 
B(k M ,z = 0) 



1 k x k z 



B 4 ,x + 



<D 



yfakj+k 2 27r(fcU 2 + l) 



B u + B 2 , z , (12) 



1 kxky 

T^B^x = Bi y + B 2 ,y, 



^k\+k\ 
1 - Anx 



Ba,x - B\ x + B 2 x . 



(13) 



(14) 



Using V ■ B2 = 0, we can derive the field B4 A from Eqs. ( 12 
14 1. Substituting the results back into Eq. ( |12) , we have the 



magnetic fields outside the superconductor B 2 z - B s z + B v z , 
with the contribution from the magnetic source 



1 



B StZ (k 2 d,z = 0) = -a 

and the contribution from the vortex 
Bw(kM,Z = 0) = (aferf+l) 

with 



k 2 d+ - ) B Uz - ik x B u - ikyBij 



On 



2^(k 2 X+l)' 



(15) 
(16) 



a (k x , k y ) = -k z [k z k 2d + ik\ - iAn X (k 2 , + k 2 )]' 1 . (17) 
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The magnetic field outside the superconductor then is given 
by 



BiJr) 



S 



d 2 kid 
2n 



Bi, z (k2d,z = 0)exp(-k 2d z+ik 2d -r 2d ). (18) 



Since the magnetic moments couple directly to the magnetic 
induction B, we use the definition that x - M x /B x . In lit- 
eratures, for examples see Refs. [T4lfT6l another definition 
X' = M x /H x was used, where H = B - 4ttM is the exter- 
nal field "seen" by the magnetic moments. The relation be- 
tween^ and^' is^-' = ^/ (l - 4nx). If x' is introduced, one 
should replace yl — 4nx in the results of the present work 
by 1 / sJl + 4nx'- Please note that the magnetic susceptibil- 
ity x = M x /B x is smaller than 1 /(4tt), i.e. x < 1 /(4tt). The 
magnetic fluctuations (M X M X ) ~ xliX ~ 4nx) diverges when 
X — > l/(47r), which indicates that the magnetic system be- 
comes unstable^. 



B. Magnetic moment perpendicular to the interface 

The calculations in this case are in parallel to those in the 
previous section. Here we skip the detailed calculations and 
only present the final results. The results can be obtained from 
Eqs. ( fT3] l and ( [To} by replacing a and k z with 



a (k x , k y j = - (k 2d + ik z ) 1 , 



1-2 



14(1-4^). 



(19) 



(20) 



IIL APPLYING TO MFM AND SSM 

We have derived the general expressions for the magnetic 
field distribution B 2z outside the superconductor in response 
to the source field B\ ,. In the following subsection, we con- 
sider the cases of MFM and SSM respectively. For MFM, the 
magnetic source is modeled as a point dipole or monopole, 
and we then calculate the force between the MFM tip and su- 
perconductors. For SSM, the source is modeled as a current 
loop and we calculate B 2 z . In both cases, the source magnetic 
field is extremely weak thus no additional vortex is induced 
by the source. 



A. Magnetic force microscopy 

In MFM, the force between the magnetic tip and the su- 
perconductor is measured as function of the distance between 
therrP2l. To calculate the force, one needs to know the mag- 
netic field distribution inside the tip. Theoretical modeling 
of the tip is challenging since the magnetic field distribution 
and shape of the tip are generally unknown. In most treat- 
ments one assumes a single cylindrical magnetic domain with 
spatially uniformly distributed moments perpendicular to the 
sample surfacepQEUlf the length of the cylinder is much larger 



than its radius, one can approximate the tip as a magnetic 
monopole. Otherwise the tip behaviors as a dipole. First we 
model the MFM tip by a point dipole along the z direction 



M, = mod (x) 6 (y) 6(z - a)z, 



(21) 



where a is the separation between the MFM tip and the su- 
perconductor. The approximation of the tip by a point dipole 
is valid when the size of the tip is much smaller than a. The 
typical size of the tip is tens of nanometer. For a » Ai, it was 
shown that the shape of the MFM tip will not affect the results 
substantially. 5 Bj then can be expressed as 



B\jy. 2 d,z = 0) = m exp(-ak2d)k2d, 

Bu(k M ,z = 0)\ . /Jt x \ 
B 1 ^ u ,z = Q)) = imeXV{ - ak7d) \k y )- 

The interaction between the tip and magnetic field is 

U(a) = - J d 3 rM z (B u + B 2 , z ) 
= ~m [B u (0,0,a) + B 2 , z (0,0,a)]. 



(22) 



(23) 



(24) 



The force then is given by F 
contribution from the source 



-d„U(a) = F s + F v with the 



'-IS 



(2ak 2c i + 1) exp(-2ak 2( j)]i2 d d k 



'<2d. 



(25) 



and the contribution from the vortex 
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exp(-ak 2d )d k 2d . (26) 



Analytical expression for the force can be obtained when 
a » Ai. In this case, only small contributes to the integra- 
tion. For the magnetic moments parallel to the interface, we 
obtain 



F \\,s 



3m o 
4Ai 



( 14 
L 

4 



2 5 V1 -4n X ) 



mo<J> 



A\n 



3A 4 , 



i^\-4n X 



•,3 



(27) 



(28) 



For the magnetic moments perpendicular to the interface, 
we have 



3m 2 n (A 4 , 



4A{ 



l5\ 



A\n 



{^1 

,4 



(29) 



(30) 



The exclusion of the magnetic flux by the superconductor 
gives rise to repulsion between the tip and superconductor, 



which is described by the first term in Eqs. (27 1 and (29 1. The 
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FIG. 2. (color online) Force between a magnetic tip and magnetic 
superconductor in the absence of vortices obtained by numerical in- 
tegration of Eq. \25\ when the magnetic moments are parallel to 
the surface. For small a/A L , the interaction becomes attractive for 
nonzero ^ as shown in the inset. At large ajX L » 1, the curves are 
described by Eq. < |27| >. 



on the vortex position relative to the tip. The position de- 
pendence of the force can be readily evaluated by replacing 
Oo <— (Do exp[i(k x x v + k y y v ]. Here r v = (x v ,y v ) is the coor- 
dinate of the vortex core. For a » Ai, the vortex-position- 
dependent force is 



F(r v = 0) - F(r v ) = - 



3m O A l 2 



(31) 



To extract the penetration depth in experiments, one mea- 
sures the force as a function of a in the absence of vortex. One 
then obtains the penetration depth by fitting to theoretical ex- 



pressions, such as Eq. (25 1. In the magnetic superconductors, 
the force depends on the orientation of the magnetic moments 
with respect to the surface. For the magnetic moments paral- 
lel to the surface, one extracts an effective penetration depth 
&lI \ ^nx, see Eq. (27 1. For the moments perpendicular 



to the surface, the bare penetration depth Ai is extracted, see 



Eq. ( 29 1. By measuring the force in two different orientations, 
one can extract both x an d ^l- 

We proceed to model the magnetic tip as a monopole and 
calculate the interaction force. The magnetic field outside the 
superconductor is given by 



force does not depends on the direction of the point dipole. 
For a magnetic superconductor, the polarization of magnetic 
moment reduces energy and causes attraction, as described by 
the second term in Eqs. p7| ) and ( |29| . When the separation a 
reduces, the attraction may be even larger than the repulsion, 



as shown by direct numerical integration of Eqs. (25 1 in Figs. 
[2]and[3] The attraction increases with^. 

The interaction between the vortex and tip depends on the 
direction of the dipole and it is attractive at large separa- 
tion a » Ai when they are parallel. To visualize a vor- 
tex, one scans the tip in experiments, and the force depends 



V B = 4-7m S(x)6(y)6(z -a), V x B = 0. (32) 
Here hq is a magnetic charge. The magnetic field distribution 



outside Z?2,j can still be calculated with Eqs. (15 1 and (18 1 



with a and k z given by Eqs. ( 10 1, ( 17 1, ( 19 1, (20i depending 



on the orientation of the magnetic moments inside the super- 
conductor. Bj in this case is given by 



BuQl2d,Z = 0) = -n exp(-ak 2d ), 



fiu-(k 2 ,i,z = 0) 
B Uv (k 2d , Z = 0) 



. exp(-ak2d) I k x 



(33) 



(34) 
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The z-component force is given by F — ho#2,z(0, 0, a). When 
the magnetic moments are parallel to the surface, we obtain 
the force due to the source Fg s and the force due to vortex 
F|| )V in the limit a » Ai 
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(35) 



(36) 



For the magnetic moments perpendicular to the surface, we 
have 
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FIG. 3. (color online) The same as Fig. |2]but for the magnetic mo- 
ments perpendicular to the surface. Results are obtained by numeri- 
cal integration of Eq. <|25 [>. 
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In the limit ^ — > and Ai — > 0, one can model the su- 
perconductor as a perfect magnetic conductor and the mag- 
netic field outside the superconductor can be obtained with 
the image method. If one models the MFM tip as a mag- 
netic dipole, the repulsion between the image dipole and tip 
is 3m^l(Aa 4 ). If the tip is treated as a monopole, the repul- 



sion force is «2/(4a 2 ). Both Eqs. (27 1 and (35 i reproduce the 
limiting results. 



B. Scanning SQUID microscopy 

In SSM, one applies external magnetic field through a field 
coil and then measures the magnetic field above the supercon- 
ductor through a pickup loopR^. We model the field coil by 
a loop with current I and radius ro. Please note that when 
a » ro, the magnetic field induced by the current loop is re- 
duced to the point dipole discussed in the previous section. 
The source magnetic field due to the current loop is then given 
by Bi = V<p with the magnetic potential 



<p(r) 



! 



d k ld ^ ld ) exp [zk 2 ,i • T2d + hdZ] ■ 



r I 



<p (k 2 </) = -j— exp (-k 2d a) J { (k 2d ro) . 
ck 2 d 



(39) 



(40) 



In this case, the magnetic field outside the superconductor due 
to the source field is 



1 \ 2nla 



Ji (kidro) exp (-k 2d a) 
(41) 



and the vortex contribution is the same as Eq. ( 16 1. Here a 
depends on the orientation of the magnetic moments. For the 
moments parallel to the surface, a is given by Eq. ( 17 1 and for 



the moments perpendicular to the surface it is given by Eq. 



( 19 1. When the magnetic moments are parallel to the surface, 



the magnetic field at the center of the pickup loop is given by 
for a » A L 



B s , = 



4ca 3 



-1 + 



3A, 



' V 1 - 4 **J 



(42) 



O 



' * 2nd 1 



1 - 



2A, 



(43) 



For the moments perpendicular to the surface, the results are 



the same as those in Eqs. (42 1 and (43 i, but without the factor 



-y/ 1 - similar to the case of MFM. 

In the case of a magnetic superconductor with isotropic 
magnetic structure as studied in Ref. [14] the magnetic field 



outside the superconductor is given by Eqs. (42i and (43 1 



The extracted penetration depth from SSM measurements is 
Ail ~~ which is larger than the bare Al- This is differ- 
ent from the effective penetration depth in magnetic supercon- 
ductors Ai JI - Anx, which is smaller than the bare Ai. 
IV. CONCLUSION 



We have calculated the magnetic fields outside a magnetic 
superconductor when a magnetic source is placed on top of 
the superconductor. For the magnetic moments parallel to the 
surface, the resulting magnetic field distribution depends on 
an effective penetration depth Ail -y/1 - Anx when the distance 
between the magnetic source and superconductor is much 
larger than Ai, while for the moments perpendicular to the 
surface, it depends on Ai. The results in the present work 
can be used to measure both the susceptibility x an d the pene- 
tration depth A L in magnetic superconductor by the magnetic 
force microscopy or scanning SQUID microscopy. This can 
be achieved by changing the orientation of the crystal. 
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